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Abstract 

Linear Programming (LP) is a powerful mathematical technique 
that can be used as a tool in Life Cycle Assessment (LCA). In the 
Inventory and Impact Assessment phases, in addition to calcu¬ 
lating the environmental impacts and burdens, it can be used 
for solving the problem of allocation in multiple-output sys¬ 
tems. In the Improvement Assessment phase, it provides a sys¬ 
tematic approach to identifying possibilities for system improve¬ 
ments by optimising the system on different environmental ob¬ 
jective functions, defined as burdens or impacts. Ultimately, if 
the environmental impacts are aggregated to a single environ¬ 
mental impact function in the Valuation phase, LP optimisation 
can identify the overall environmental optimum of the system. 
However, the aggregation of impacts is not necessary: the sys¬ 
tem can be optimised on different environmental burdens or 
impacts simultaneously by using Multiobjective LP. As a result, 
a range of environmental optima is found offering a number of 
alternative options for system improvements and enabling the 
choice of the Best Practicable Environmental Option (BPEO). 
If, in addition, economic and social criteria are introduced in 
the model, LP can be used to identify the best compromise solu¬ 
tion in a system with conflicting objectives. This approach is 
illustrated by a real case study of the borate products system. 


Keywords: Allocation, multiobjective linear programming; Best 
Practical Environmental Option (BPEO); Boron; BPEO; dual 
values, LCA; LCA; Life Cycle Assessment, LCA; marginal analy¬ 
sis, LCA; operations research, LCA; Pareto Analysis; system 
optimisation, LCA 


1 Introduction 

One of the main applications of Life Cycle Assessment (LCA) 
is comparison of possible modifications to an existing prod¬ 
uct or process with the aim of improving its environmental 
performance. In most cases there are many possibilities for 
improvements and the choice of the best alternative is not 
obvious. Therefore, to identify the optimum options for sys¬ 
tem improvements, it is necessary to use a suitable optimi¬ 
sation technique. This paper illustrates the potential of ap¬ 
plying one such technique - Linear Programming - to 


optimise the environmental performance of a product sys¬ 
tem as a part of the Improvement Assessment phase of LCA. 
In addition, it demonstrates the value of Linear Program¬ 
ming (LP) in solving the problem of allocation in multiple- 
function systems as a part of both the Inventory and Impact 
assessment phases, and shows how the approach can be used 
to identify the Best Practicable Environmental Option for a 
process or product system. The use of the approach is illus¬ 
trated by a specific example of a multi-product system which 
produces different borate products from boron ores. 

2 Linear Programming and Life Cycle Assessment 

For present purposes, a model means a set of mathematical 
relationships which describe the operation of the unit proc¬ 
esses forming a product system. LCA is based on linear ho¬ 
mogeneous and unconstrained models of human economic 
activities and of their effect on the environment; i.e. envi¬ 
ronmental burdens and their impacts are assumed to be di¬ 
rectly proportional to the number of functional units pro¬ 
duced (e.g. Heijungs, 1992; Huppes and Schneider, 1994). 
The commonest and simplest form of linear model, usually 
known as an input/output model, takes the form of a set of 
homogeneous linear equations relating the outputs from and 
inputs to an economic system. In LCA, this type of model is 
extended to relate the environmental burdens associated with 
a product system to its economic outputs (e.g. Heijungs, 
1997). However, the simple input/output approach is not 
able to account for the internal structure of the system and 
it cannot help to improve the efficiency of industrial opera¬ 
tions. Simple input/output modelling is also not applicable 
to the analysis of a product system whose operation is con¬ 
strained, for example, by the capacity of existing unit proc¬ 
esses, or by the availability of material or energy inputs. 
Since in reality, almost all systems are subject to a number 
of constraints, a more powerful approach to system analysis 
and optimisation in LCA is necessary. 

The approach proposed here uses Linear Programming (LP) 
to model the behaviour of a linear system subject to con¬ 
straints. The same method can be used for non-linear sys- 
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terns that can he approximated as linear. The general math¬ 
ematical relationships are given in the text, with a more de¬ 
tailed account of LP in the Appendix. Other examples of the 
applications of LP to LCA, in addition to the specific exam¬ 
ple shown here, are given by Azapagic (1996) and Azapagic 
and Clipt (1994,1995a, 1995b, 1998a, 1998b, 1998c). The 
concepts behind the use of LP in LCA are shown schem¬ 
atically in Figure 1. 

The main characteristic of this kind of modelling is that it is 
based on physical and technical relationships between the 
inputs and outputs and environmental interventions of the 
system. Therefore, a LP model describes the underlying physi¬ 
cal causation in the system and thus lends itself naturally to 
solving allocation in multiple-function systems according to 
the procedure recommended by ISO 14041 (1997). Moreo¬ 
ver, because LP modelling describes complex interactions 
between different parts of the system, it can describe changes 
in the operating state of the system and associated environ¬ 
mental interventions, resulting from changes in material or 
process properties. This approach therefore reveals how en¬ 
vironmental burdens and impacts - and their allocation be¬ 
tween different functions - change as the operation of the 
system is changed. These features are particularly useful in 
the Inventory and Impact Assessment phases, as discussed 
further in the paper. 

The LP approach is also valuable in the Improvement As¬ 
sessment phase for quantifying the trade-offs between dif¬ 
ferent environmental burdens and impacts. In effect, the LP 
model identifies a range of optimum solutions for improved 
environmental performance which can be achieved through 
various modifications to the product system, including 
changes in economic performance. In this way, by combin¬ 
ing both environmental and socio-economic criteria, this 
approach enables identification of the Best Practicable Envi¬ 
ronmental Option (BPEO) not entailing excessive costs. 



Fig. 1: Linear Programming as a tool in Life Cycle Assessment 


2.1 Modelling in the Linear Programming Format 
A conventional Linear Programming model of an economic 


system has the form: 

Maximise (or Minimise) F = 

if.*. 

i=i 

0) 

subject to 

l 

E a c.i X i - e c 

i=l 

c= 1,2,...,C 

(2) 

and 

_X 

IV 

O 

i = 1,2,...,1 

(3) 


where eqn. (1) represents an objective function, usually a 
measure of economic performance (e.g. profit or cost) and 
eqns. (2)-(3) are linear constraints in the system, describing 
material and energy balance relationships, productive ca¬ 
pacity, raw material availabilities, quality requirements, 
market demand and so on. The constraints can be defined 
as equalities or inequalities. The LP example in the Appen¬ 
dix shows how mass and energy balances lead to equality 
constraints, whereas raw material availability and other con¬ 
straints can lead to inequalities. The variables x. represent 
quantitative measures of material and energy flows includ¬ 
ing inputs, flows within the economic system, and outputs. 
In LP these variables are termed "activities", or more pre¬ 
cisely, "activity levels". The coefficients a t . are the factors of 
proportionality between activity (process) levels and inputs 
and outputs from the system. The right hand side coeffi¬ 
cients, e , represent the limitations on the constraints. 

The optimum point of a LP model is defined by the "active" 
constraints, i.e. by the constraints that are satisfied as 
equalities at the solution of the LP model. The other, non¬ 
active constraints do not influence the solution of the sys¬ 
tem and usually have "slack" or unused resources associ¬ 
ated with it (—> Appendix, p. 314). At the solution, each 
active constraint has a dual or marginal value, which shows 
the change in the objective function with the change in the 
right hand side coefficient of that constraint. The discussion 
in the following section shows that the dual values represent 
the allocation coefficients which describe the physical cau¬ 
sality in the system. Therefore, LP is useful in the Inventory 
and Impact Assessment phases for solving the allocation 
problem that arises in multiple-function systems. 

The types of questions that conventional LP helps answer 
are, for instance, related to finding the optimum operating 
point in the system that maximises the profit (or minimises 
the costs) and uses the optimum amount of resources, sub¬ 
ject to the constraints. In the context of LCA, the general LP 
model has the same form; however, the constraints (2) now 
encompass all activities from extraction of primary materi¬ 
als from earth through processing to final disposal. In addi¬ 
tion, the functional output or outputs are also treated as 
activities. Furthermore, the objective functions are now de¬ 
fined by the environmental burdens if the analysis is at the 
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Inventory level, rather than an economic objective, as repre¬ 
sented by Azapagic (1996): 

Minimise B = V be . x. (4) 

i Z-i J. > ' 

i = i 

where be. is burden j from process or activity x.. The objec¬ 
tive functions can also be defined as the environmental im¬ 
pacts: 

J 

Minimise E k =^ec k jBj (5) 

j=1 

where ec k . represents the relative contribution of burden B ; 
to impact E t , as defined by the "problem oriented" approach 
to Impact Assessment (Hkijungs et ah, 1992). The LP exam¬ 
ple in the Appendix also show how environmental objec¬ 
tives can be formulated. 

Depending on the goal of the study, the system can be 
optimised on one or a number of environmental and eco¬ 
nomic objective functions, to identify the optimum solutions 
for system improvements. This application of LP in Improve¬ 
ment Assessment is explained below. 

2.2 Linear Programming in Inventory and Impact 
Assessment 

The aim of the Inventory and Impact Assessment phases is 
to quantify total environmental burdens and impacts and 
not to optimise performance of the system. Therefore, the 
optimisation is not performed at this level of analysis; the 
LP model is solved to calculate the burdens and impacts for 
possible operating states of the system, which are normally 
defined by the optimum economic performance. For the type 
of systems that most LCA studies consider, i.e. unconstrained 
and homogeneous, this step is equivalent to the conventional 
LCA procedure of calculating the burdens and impacts. 
However, unlike the simple input/output analysis in the con¬ 
ventional approaches, formulation of the system as a LP 
model also provides a rigorous basis for allocation of the 
burdens and impacts in multiple-function systems by physi¬ 
cal causality. This approach has already been proposed by 
Azapagic and Clift (1994, 1995a, 1995b, 1998b) and 
Azapagic (1996) and is explained in more detail below. 

2.2.1 Allocation in multiple-function systems 

The problem of allocation in multiple-function systems is to 
find a procedure which would assign to each of the func¬ 
tions of the system only those environmental burdens and 
impacts for which it is responsible. There are three types of 
multiple-function systems where allocation can be relevant: 
multiple-input systems (waste treatment processes), multi¬ 


ple-output systems (co-production), and multiple-use sys¬ 
tems (open-loop recycling). Because the method for solving 
the allocation problem may depend on the goal of the study, 
different approaches to allocation are possible: allocation 
may be avoided by enlarging or disaggregating the system, 
or it can be solved by applying an appropriate causation 
principle (ISO, 1997). Avoiding allocation by system enlarge¬ 
ment is an appealing way to deal with this problem; how¬ 
ever, in many cases the systems considered will become more 
complicated or uncertain because of the additional data 
needed. Avoiding the allocation problem by disaggregation 
may be possible, but only if detailed data about the system 
are available and if the system can be broken down to elimi¬ 
nate all processes common to the functional outputs. 

In all other cases some kind of allocation will still be neces¬ 
sary. The recommended approach (ISO, 1997) is to base al¬ 
location on the physical 1 causal relationships between bur¬ 
dens and functional outputs, which in turn requires a model 
of the system behaviour. Therefore, allocation on an arbi¬ 
trary basis, such as mass or molar flow, must be avoided, 
unless the modelling of causalities shows such a basis to be 
appropriate. However, if it is not possible to apply the cau¬ 
sation principle, for example if the functional outputs can¬ 
not be varied independently (see below), then economic re¬ 
lationships should be used as the only other relevant basis 
for allocation (Hupphs and Frischknkcht, 1995; ISO, 1997; 
Clift et ah, 1998). 

When the causal relationships are represented by a model 
which describes the real behaviour of the product system, 
the model can be used to allocate burdens between different 
functions by exploring how the burdens change when the 
quantity of one function is changed with the quantities of all 
the other functions kept constant (e.g. Azapagic and Ci.iit, 
1998b). An example of this kind of problem in the chemical 
industry is naphtha cracking, where the outputs can be var¬ 
ied independently by changing cracking conditions (Azapagic 
and Clift, 1995b). 

The type of changes considered can be marginal, incremen¬ 
tal or average which, in turn, depends on the goal and scope 
of the study and the questions to be answered by LCA (Clift 
et al., 1998; Azapagic and Clift, 1998b). Marginal changes 
are relevant when the performance of a specific system is 
analysed to determine the effect of infinitesimal changes. 
Because the changes are infinitesimal, this amounts to ana¬ 
lysing the present state of the system, which most LCA stud¬ 
ies are initially concerned with. This kind of analysis is rel¬ 
evant if the goal of the study is comparison of the products 
from a multiple-function system or comparison of one of 
the products from this system with a product from another 
system. Marginal changes can always be described by a lin¬ 
ear model, even if the model is not homogenous (Clift and 


1 The term "physical relationships" has a broader meaning in this context 
and includes physical, chemical, biological and technical relationships. 
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Azapagic, 1995; Clift et al., 1998). Therefore, LP model¬ 
ling can always be used to describe marginal changes. As 
discussed in the following section, the marginal values cal¬ 
culated at the solution of a LP model represent the marginal 
allocation coefficients which reflect the physical causal rela¬ 
tionships in the system. 

However, if the goal of the study is to consider incremental 
or average changes 2 to a new state of the system, for instance 
to describe product or process development, then a linear 
model may not be applicable (Clift and Azapagic, 1995; 
Clift et al., 1998). Where the system behaviour can be 
linearised about the alternative operating states being con¬ 
sidered, it may still be possible to use a LP model. This is 
normally the case for incremental changes, but is less likely 
to be applicable for average changes. 

However, in some systems, the ratio between two or more 
functional units and their parameters in the system may be 
fixed; examples of this arise in the chemical industry, where 
the ratio of sodium-hydroxide (NaOH) to chlorine (Cl,) pro¬ 
duced by electrolysing brine is constant. In this case, alloca¬ 
tion by physical causation cannot be used and socio-eco¬ 
nomic relationships, as the only other relevant choice, should 
be used instead. The argument for this is that economic re¬ 
lationships reflect socio-economic demand which causes the 
multiple-functional units to exist at all (Clift et al., 1998). 
These relationships are usually translated into some economic 
measure, such as net value at the point in the system where 
the product streams divide (Huppes and Frischknecht, 1995). 


or in a simplified notation: 

C 

dF=XMec <8> 

C=I 


The partial derivative: 



(9) 


represents the dual or marginal value which is calculated 
automatically at the solution of a LP model. For a marginal 
change, the derivative (9) will remain constant, i.e. the state 
of the system will not change, so that after integration, eqn. 
(8) becomes: 

F=tVc < 10 > 

C=1 


so that, in these terms, the model is both linear and homog¬ 
eneous. For the objective function defined as environmental 
burden B, then by analogy with eqn. (10), the total value of 
that burden is related to the dual values by: 



2.2.2 Linear programming and allocation 

The solution of the LP model defines for each constraint, as 
given by eqn. (2), a dual or marginal value which shows the 
effect on the objective function of a marginal change in the 
right-hand side coefficient, e , of the constraint with all the 
other constraints unchanged. This analysis is applicable only 
where the coefficients e |V ..e c can in principle be subject to 
independent marginal changes. 

The value of the objective function, as defined by eqn. (1), 
at the solution of the model can then also be written as a 
function of the right-hand side coefficients: 

F = /[e,,e 2 ,...,e c ] (6) 


where l c is the marginal or dual value of the cth constraint 
for burden j, i.e. it represents the change in the total burden 
Bj with the change in coefficient e .: 


J.c 


3B: 

de c 


. e c 


( 12 ) 


Because the environmental impacts are linear homogenous 
functions of the burdens, they can be formulated in the same 
way. In general, eqn. (11) can be written as: 


c 

^k — 


c=l 


(13) 


In the case of marginal changes in these coefficients, the cor¬ 
responding change in the objective function is equal to: 


c=l ac c 


(7) 


2 Incremental changes are small but finite changes which change the operat¬ 
ing state of the system. Average changes are related to more substantial 
changes, for instance eliminating a functional output completely. 


where: 


M%c = 


/3E k 

3e c 


'e. . .. e c 


(14) 


The LP model can be formulated so that the quantities of 
functional outputs are represented by coefficients e c , amongst 
other coefficients describing constraints of the system. If at 
the solution of the LP model only those constraints that are 
related to the functional outputs are active, then the dual 
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values represent the marginal allocation coefficients, which 
relate changes in a burden or impact to the infinitesimal 
changes in one output while all other functional outputs are 
held constant; i.e. 1 and m ki; represent the allocated burden 
or impact, respectively. In LCA terms, such a system is un¬ 
constrained so that the burdens and impacts are allocated 
fully between the outputs (Azapagic, 1996; Clift et al., 1998; 
Azapagic and Clift, 1998b). Naturally, as noted above, the 
dual or marginal analysis is only applicable when all the 
coefficients e c can be varied independently. Thus, as expected, 
the marginal allocation approach cannot be used for out¬ 
puts which are produced in fixed proportions. However, 
whenever constraints which describe capacity or market 
demand are active at the solution, the burdens and impacts 
will be allocated in part to these constraints so that they are 
not fully allocated to the functional outputs (—» Appendix). 
This is consistent with the physical causation principle, as 
the changes in the system behaviour in this case are not de¬ 
termined by the functional outputs only, but by other prop¬ 
erties of the system as well. 

An illustration of the marginal allocation approach for a 
multiple-function system producing five boron products is 
presented in Section 3. 

2.3 Linear Programming in Improvement Assessment 

The main objective of the Improvement Assessment phase is 
to identify opportunities for improving the environmental per¬ 
formance of the system. To achieve this, the system has to be 
optimised on the environmental objective functions, B. or E k . 
There are different ways to approach the optimisation prob¬ 
lem, depending on the goal of the study. For instance, the en¬ 
vironmental impacts could be aggregated into a single impact 
function in the Valuation phase and the system optimised on 
this objective to give the overall environmental optimum of 
the system. However, we argue that, until the methodology is 
developed further. Valuation should be avoided and the sys¬ 
tem optimised on the set of environmental burdens or impacts 
instead (Azapagic, 1996). This is possible by using 
Multiobjective LP (MOLP) which enables simultaneous 
optimisation on a number of objective functions, resulting in 
a range of environmental optima of the system. These optima 
define the multidimensional non-inferior or Pareto surface, 
which is optimal in the sense that none of the objective func¬ 
tions can be improved without worsening one or more of the 
other objective functions. Therefore, trade-offs between ob¬ 
jective functions are necessary in order to select the best com¬ 
promise solution. For example, if two objective functions, de¬ 
fined as environmental impacts, are optimised simultaneously, 
the resulting Pareto optimum does not necessarily mean that 
these impacts are at their minima obtained when the system is 
optimised on each of them separately. The Pareto optimum, 
however, does mean that the set of best possible options has 
been identified for a system in which both impacts should be 
reduced. The value of MOLP in LCA, therefore, lies in offer¬ 
ing a range of alternative solutions; they are all optimal in the 


Pareto sense, but the choice of the best one will depend on 
preferences and constraints imposed on decision-makers. 

However, decisions are rarely made on the basis of environ¬ 
mental LCA only and other considerations, particularly eco¬ 
nomic and social, are usually involved. If, in addition to the 
environmental objectives, multiobjective optimisation also 
includes economic and social objectives, acceptable solutions 
which represent the compromise between conflicting objec¬ 
tives can be identified. This approach could be of particular 
importance for the process industries, which face the prob¬ 
lem of having to keep production costs down while at the 
same time complying with environmental legislation and 
responding adequately to increased public awareness of en¬ 
vironmental problems. For example, it shows how far from 
its economic optimum the system must be operated to achieve 
certain improvements in environmental performance. 

The application of MOLP to the Improvement Assessment 
phase is now illustrated by a real case study of a multiple- 
function (co-product) system producing five borate products. 

3 Case Study 

The co-product system considered here produces the follow¬ 
ing products: 5 and 10 mol borate, boric acid (BA), anhy¬ 
drous borax (ABA) and anhydrous boric acid (ABA). A sim¬ 
plified flow diagram of the process is shown in Figure 2. 
Further details of the system are given elsewhere (Azapagic, 
1996; Azapagic and Clift, 1998d). Boron minerals - borax 
and kernite - are extracted in the mine, crushed and trans¬ 
ported to the plant. 5 and 10 mol borates are produced by 
dissolving borax and kernite in water; Na-borates are then 
separated from insolubles, crystallised and dried to produce 
powder products. Boric acid is produced in a separate plant, 
by reacting kernite ore with sulphuric acid. The rest of the 
process is similar to the 5 and 10 mol production. Anhy¬ 
drous borax and anhydrous boric acid are made in high- 
temperature furnaces from 5 mol borate and BA, respec¬ 
tively. All products are then either packed or shipped in bulk. 
Electric energy and the steam for the system are provided by 
the on-site natural gas cogeneration facility, which meets all of 
the electricity and most of the steam demand. Any additional 
steam is provided by a steam plant which is also fired by natu¬ 
ral gas. The overburden from the mine and the gangue from 
the process are stocked in piles; the waste water from the re¬ 
finery is discharged into self-contained ponds. All activities, 
from extraction of raw materials to the production of the bo¬ 
ron products and materials used, are included in the system. 
However, the use and disposal phases of the products are not 
considered in this study ("cradle-to-gate" approach). 

The goal of the study is to evaluate the environmental per¬ 
formance of the system and to identify the opportunities for 
environmental improvements. Since the study is intended for 
internal use, the functional unit is taken to be total yearly 
production of the boron products, i.e. 1062000 t/yr. 
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Fig. 2: LCA flow diagram of rhe boron products system 


3.1 Allocation in the boron products system 

The LP model of the boron system is described by the con¬ 
straints defined by material balances, output of products, 
supply of primary and raw materials, capacity and heat re¬ 
quirements as given by: 


I 

Mass balance: ^a c jXj = 0, i=l,2,...,l (15) 

i=l 


by the supply S m , activities or processes x i are subject to the 
capacity limit Cand the heat requirement H y is constrained 
by the heat availability Q,. Since the discussion in this chap¬ 
ter is related to the functional unit defined as the operation of 
the system for one year, the product demand D t| is taken to be 
equal to the total output of each product for one year. The 
objective functions of the system are environmental burdens 
and impacts, as defined by eqns. (4) and (5). A large scale LP 
software, XPRESSMP (Dash Associates, 1993), has been used 
to solve the system model and to calculate marginal values. 

The functional outputs from the system can be changed in¬ 
dependently and therefore allocation based on physical cau¬ 
sation is appropriate. Since this study is concerned with a 
specific system and marginal changes around its operation, 
the burdens and impacts are allocated through marginal 
values of the LP model, found at the solution of the model. 
As noted above, for the purposes of allocation, rhe model is 
not optimised at this stage; it is solved to calculate the total 
burdens and impacts and the corresponding marginal val¬ 
ues at the operating point of interest. The solution is found 
at the intersection of active constraints, which are in this 
case defined by the quantities of the product outputs (see 
eqn. (16)). Therefore, these constraints have non-zero mar¬ 
ginal values, which are equivalent to rhe marginal alloca¬ 
tion coefficients. Other constraints are non-active and their 
marginal values are zero; therefore, the burdens are fully 
allocated to the products. If the operating state of rhe sys¬ 
tem changed so that other constraints determined the op¬ 
eration of the system, then they would also be active at the 
solution and their marginal values would be non-zero; hence 
the burdens would be allocated to these constraints as well 
as to the products. 

As an example, the results of the marginal allocation of the 
total CO, emissions of 295860 t/yr among different prod¬ 
ucts are shown in Table 1 and Figure 3. For instance, the 
marginal value, i.e. the emission of CO, allocated to 10 mol 
on the marginal basis, is equal to 0.18 t/t; i.e. if the output 
of this product increases by 1 t, the total CO, emissions 
increase by 0.18 t. The CO, emissions allocated to other 
products range from 0.24 t/t for 5 mol, through 0.41 t/t for 
BA to 0.98 t/t and 1.94 t/t for AB and ABA, respectively. 


Products output: 

P q <D q ,q=l,2,...,Q 

(16) 

Primary and raw 
material supply: 

R nl <S m ,m=l,2,...,M 

(17) 

Productive capacity: 

X; <Ci,i=l,2,...,l 

(18) 

Heat requirements: 

H z < Q z , z=l,2,...,Z 

(19) 


where production P q is limited by the product demand D , 
primary and raw materials consumption R m is determined 


Table 1: CO, emissions allocated to the boron products 


Product 

Functional 

Output* 

(t/yr) 

Marginal 

allocation 

(t/t) 

Total 
allocated 
CO z (t/yr) 

10 mol 

81000 

0.18 

14580 

5 mol 

810000 

0.24 

194400 

BA 

150000 

0.41 

61500 

AB 

16000 

0.98 

15680 

ABA 

5000 

1.94 

9700 

TOTAL 

1062000 


295860 


'Production figures shown in this table are not based on actual 
production or sales figures in any particular year, but represent a 
possible production scenario. 
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Let us now compare the results of the marginal allocation 
approach with the arbitrary allocation methods most com¬ 
monly used in co-product systems: mass and market value 
bases (—> Fig. 3). To illustrate the importance of system dis¬ 
aggregation, allocation on the mass basis for both aggre¬ 
gated and disaggregated systems is considered. System dis¬ 
aggregation was possible in this particular case study be¬ 
cause the detailed data for the system were available. If allo¬ 
cation on the mass basis is done without system disaggrega¬ 
tion, the total burdens are allocated among different prod¬ 
ucts in proportion to the mass of their outputs. In that case 
all products in the boron system would have the same allo¬ 
cation factor of 0.28 t/t (total CO, emissions of 295860 t/yr 
divided by the total production of 1062000 t/yr). For the 
output of 10 mol of 81000 t/yr, for instance, the CO, emis¬ 
sions allocated to this product would be equal to 22680 t/yr. 
The marginal approach allocates 14580 1 CO,/yr to the same 
product, which represents a total difference of 36%. For 
ABA the difference is much higher and is equal to 86%. The 
reason for this is that allocation by mass without system 
disaggregation implicitly assumes that all processes for pro¬ 
duction of co-products are the same, which is obviously not 
true. These results illustrate the point that allocation done 
in this arbitrary way does not reflect physical causality and, 
therefore, cannot allocate the burdens realistically in com¬ 
plex industrial multiple-function systems. A similar con¬ 
clusion emerges if allocation is based on the B,0, content 
in different products. Although it seems that the B,0, allo¬ 
cation coefficients increase as the degree uf processing of 
different products increases to give a higher boron content 
(—> Fig. 3), this increase is actually related linearly to the 
content of B,0, in the products and not to the burdens asso¬ 
ciated with their processing. 


that in case of the CO, emissions, allocation based on a physi¬ 
cal quantity, i.e. mass, is appropriate. Thus, it can be correct 
to allocate the burdens on the basis of a physical quantity, 
provided that this basis emerges from analysis of causation 
in the system. 

Allocation on the basis of market value’does not give cor¬ 
rect results in this case. Although the burdens allocated on 
the marginal and market value bases are quite similar for 
lOMol, 5MoI and BA, the difference is much larger for AB 
and ABA (51% and 38%, respectively). This implies that 
the external costs of the environmental burdens are not pro¬ 
portional to the current economic values of these products. 
Therefore, allocation by financial value can give misleading 
results and should not be used in systems where physical 
causality exists. 

Allocation by physical causality and, in particular, the use 
of LP for marginal allocation are discussed in more detail in 
Azapagic and Clut (1998b, 1998c). 

3.2 Improvements in the boron system 

In addition to evaluating the environmental performance, 
the other goal of performing a LCA study of the boron prod¬ 
ucts system was to identify a range of possibilities for mini¬ 
mising total environmental impacts from the system, while 
maximising production subject to product demand and keep¬ 
ing the production costs at a minimum. The information 
obtained would then serve as a basis for effecting improve¬ 
ments in the system. 


However, the situation is quite different if allocation on the 

mass basis is done after the system has been disaggregated _ 

to take into account differences in the processes for produc- * The market value is here raken to be gross selling price of the boron pend¬ 
ing different products. In that case allocation on marginal ucts, respectively: lOMol: $239.5; 5Mol: $255.4; BA: $527.4; AB: $612.5; 

and mass bases gives the same results (-» Fig. 3). This means ABA: 5 1535.2 



■ Marginal basis 

ESS Mass basts (aggregated data) 

i 1 B 2 O 3 Content 

E3 Mass basis (disaggregated data) 

K Market value basts 


Fig. 3: Comparison of different allocation methods for CO, emissions 
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The objective functions of the LP model, therefore, include 
environmental impacts, costs, and total production. To en¬ 
able the approach to be illustrated graphically, the system 
is here optimised on three objectives only: Costs (F), Glo¬ 
bal Warming Potential (GWP) and Production (P) as de¬ 
fined by: 

I 

Minimise F — X fjXj (20) 

i=I 

J 

Minimise GWP = Y ec-B, (21) 

i—t j j 
j=i 

Q 

Maximise P =: ^ (22) 

q = i 

subject to the constraints defined by eqns. (15)-(19). Obvi¬ 
ously, eqns. (20) and (21) are equivalent to eqns. (1) and 
(5) and they are expressed in £/yr and t/yr, respectively. Equa¬ 
tion (22) is defined by the total output of the five boron 
products in t/yr (—> Table 1). Therefore, the model is simi¬ 
lar to that used for allocation of environmental burdens, 
except that the product output constraints (16) are now 
different: instead of being defined by current operations, 
they are determined by the market demand projected on 
the basis of the trends in the previous years. In addition, 
several alternatives to the way the system is operated at 
present have been considered to identify the Best Practica¬ 
ble Environmental Options (BPEO). These alternatives in¬ 
clude transport of the ore in the mine by a conveyor instead 
of trucks, and steam generation in the Cogeneration instead 
of the Steam plant. Depending on the position on the non¬ 


inferior surface, different alternatives are chosen in the 
optimisation and they represent the BPEO for that particu¬ 
lar operating state. As the position on the surface changes, 
so does the BPEO. This is illustrated in Figure 4. 

As a result of multiobjective optimisation on functions (20)- 
(22), the non-inferior surface ABCD with a range of opti¬ 
mum solutions is generated. The results shown in Figure 4 
are normalised to the optimum values of each objective, i.e. 
F*, GWP* and P*, obtained in single-objective optimisation, 
when the other two functions are ignored. Point A in Fig¬ 
ure 4 represents the minimum of the cost objective func¬ 
tion; however the production is at the minimum and GWP 
is 31% above its optimum value (G*). For instance, the 
BPEO at this point includes transport in the mine by the 
trucks and steam produced in the Steam plant. If a further 
decrease in the GWP objective by one tonne is required, 
that would result in a cost increase of £95 for constant prod¬ 
uct output; i.e. the marginal cost for reducing global warm¬ 
ing potential from the system in this case is £95. Similarly, 
if the production were to increase by one tonne, the result¬ 
ing increase in the costs would be £34 if GWP were kept 
constant. 

By moving from point A along the non-inferior curve for 
constant GWP, both costs and production increase, to reach 
their maximum feasible values at point B. Here, the Cost 
function is 4% above its optimum value. If Production is 
increased by one tonne, £300 of the Costs objective have to 
be given up. Similarly, one tonne change in the GWP is as¬ 
sociated with a cost change of £100; i.e. the marginal cost 
of reducing GWP changes as a result of the change in oper¬ 
ating conditions. At this solution, the BPEO is defined by 
steam generation mainly in the Steam plant and the pre¬ 
ferred transportation means in the mine are the conveyors. 
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Furthermore, if for instance the system were to be operated 
at point C, GWP would be 3.3% above its optimum value 
obtained in the single-objective optimisation. The produc¬ 
tion would be at the minimum, and the costs would increase 
by 14%. The effect of GWP and Production on Costs is 
similar to that found at point A: an improvement in GWP of 
one tonne would worsen the Costs objective by £100, while 
a tonne increase in P would result in £36 increase in costs. 
Here, 93% of the steam is generated by the Cogeneration 
plant and the rest is produced in the Steam plant. The con¬ 
veyors still remain the best transport option in the mine. 

However, if for example, point D were to be chosen as the 
best compromise solution, then for the same value of GWP 
as at point C, the production would reach the maximum; 
however, costs would have to increase by 17%. It can be 
noticed here that both GWP and Production exhibit similar 
effect on the Costs: a decrease in GWP by one tonne in¬ 
creases Costs by £3600. If the production is increased by 
one tonne, the costs increase by £3500. At this solution, the 
best practicable environmental option is defined by truck 
transport in the mine and steam production in the Cogenera¬ 
tion plant. It may be noted that, for instance, the marginal 
costs of reducing GWP at point D are by a factor of 100 
higher than at point C due to the different environmental 
options chosen at these two points. 

These results demonstrate how optimum solutions, and there¬ 
fore, BPEO, change with the operating state of the system. 
The same analysis can be carried out for other points on the 
non-inferior surface which are all optimal in the Pareto sense. 
By trading-off the values of different objectives at these 
points, decision-makers can select any solution on the sur¬ 
face, depending on how much of one objective they are pre¬ 
pared to give up in order to gain in another. If all objectives 
are considered to be of equal importance, then one of the 
possible ways to choose the best compromise solution is to 
identify operation by which all objectives differ from their 
optima by the same percentage. However, should some ob¬ 
jectives be considered more important than the others, any 
other solution on the non-inferior curve can be chosen as 
the best compromise. The value of multiobjective 
optimisation in the context of LCA, therefore, lies in offer¬ 
ing a range of choices for environmental and economic im¬ 
provements of the system and so enabling preferences to be 
identified after analysing all the trade-offs among objectives. 

4 Conclusions 

Linear Programming (LP) is a powerful mathematical tech¬ 
nique that can be successfully combined with Life Cycle 
Assessment (LCA). Marginal analysis, which is an integral 
part of LP, can be used to solve the problem of allocation in 
multiple-function systems. LP is also useful in identifying 
opportunities for environmental improvements in a product 
system. Since improvements cannot be carried out on the 
basis of environmental LCA only, LP can also be used to 


quantify the compromise between environmental and eco¬ 
nomic performance by using multiobjective optimisation. The 
advantage of multiobjective optimisation in environmental 
system management in the context of LCA lies in providing 
a set of alternative options for system improvements rather 
than a single optimum solution, and enabling the choice of 
the Best Practicable Environmental Options (BPEO) by re¬ 
vealing the relationships between environmental improve¬ 
ments and costs. 


Nomenclature 

a c,i 

Input/output or proportionality coefficients in 


constraints 

bc j.i 

Environmental burden coefficients 

1 

Environmental burdens 

C 

Total number of constraints in LP model 

q 

Capacity of activity i 

°q 

Demand for product q 

e c 

Right-hand side coefficients in constraints 

ec k,j 

Environmental impact coefficients 

E k 

Environmental impacts 

F 

(Economic) Objective function (profit or costs) 

F* 

Optimum value of the cost function obtained in 
single objective optimisation 

h 

Coefficients in economic objective function 

GWP 

Global Warming Potential objective function 

CWP‘ 

Optimum value of the GWP objective function 
obtained in single objective optimisation 

H z 

Heat requirement 

I 

Total number of activities 

J 

Total number of burdens 

K 

Total number of impacts 

M 

Total number of primary and raw materials 

P 

Total production objective function defined by the 
products output Pq 

P* 

Optimum value of the P objective function obtained 
in the single objective optimisation 

p q 

Product outputs (functional units) 

Rm 

Primary and raw materials 

Q 

Total number of products from the boron system 

Qw 

Heat lost or dissipated from process 

Q z 

Heat supply 

s m 

Supply of primary and raw materials 

x i 

Activity level (operation level of process) 

H* 

Marginal values in LP model and marginal 
allocation coefficients for burden j and constraint c 

Fk,c 

Marginal values in LP model and marginal 
allocation coefficients for impact k and constraint c 
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Appendix - An Example of a LP Model 

For readers less familiar with the mathematical concepts of 
LP, a simplified example showing how to formulate and solve 
a LP model in the context of LCA is presented here. 

Consider a multiple-function system which produces two prod¬ 
ucts and uses five alternative feedstock materials (—» Fig. A. 1 ). 
The inputs into the system are denoted by x, to x 5 and the 
outputs by x ( and x 7 . In the LP terms, the variables x, to x 7 



► * 


*• x- 


Fig. A.l: Multiple-function system with two products 


are known as activities, or more precisely, activity levels. 
For the purposes of this example, it is assumed that the sys¬ 
tem produces only one environmental burden, B, and con¬ 
sumes a certain amount of energy, Q. The system produces 
1000 kg of each product and they are taken to be the func¬ 
tional units of the system. 

The raw materials are in two categories: x, and x, from the 
first and the other three feeds from the second category. The 
amounts of these two feed categories are limited and it is 
not possible to process more than 1500 kg of the first cat¬ 
egory and more than 1200 kg of the second. The costs (per 
tonne) of the raw materials are: 

X, x, X, x 4 X, 

Cost (£/t) 110 120 130 110 115 

The product x 6 sells for £150 per tonne and x 7 for £120/t. 
The unit energy used in the system costs 0.3 pence/kg of 
each product. 

In conventional LP modelling, the most common type of 
questions asked is: How should the manufacturer make these 
products in order to maximise their profit? LCA, though, is 
not so much concerned with the economics of the system as 
with the question of how to minimise the environmental 
impacts. The rest of this example demonstrates how LP can 
provide answers to both of these questions. However, the 
model building is demonstrated first. 
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Model building 

A number of constraints can be defined according to the 
specifications given above. First of all, the mass and energy 
balance is carried out, to define the mass and energy bal¬ 
ance constraints. This step is equivalent to what is normally 
done in the Inventory Analysis phase. Therefore, we have: 

Input = Output 

Mass balance: X, + X, + X 3 + X 4 + X 5 = X 6 + X 7 + B 

In the case where emissions are small compared with the 
process flows, as is the case in this example, B can be ne¬ 
glected. In the LP terms, the above equation is normally 
written in the following way: 

Mass balance constraint: 

x 1 + x 2 +x 3 + x 4 + x 5 —X 6 -X 7 =0 (A.l) 

The mass balance constraints are therefore always defined 
as equalities. 

The energy balance around the system in its most basic form 
is defined by the following equation: 

Energy balance: Q z = H product - H feed +Q W 

where H pr(lduct and H (mJ are the enthalpies of the product and 
the feed, respectively, expressed in kj; Q y is the enthalpy used 
in the system and Q w is the enthalpy lost. As all constraints 
and objective functions in LP models are expressed per unit 
level of the activity, the enthalpies in the above equation can 
be substituted by the specific enthalpies, h,,h„...,h 7 of the 
feed the and products. They are shown below: 

x, x, x, x, x, x„ x 7 

h (kj/kg) 1200 1300 1100 1200 1300 1200 1400 

Thus, the energy balance constraint is equal to: 

Energy balance constraint: 

1200x 6 + 1400x 7 — 1200x, - 1300x 2 - 
1 lOOx 3 - 1200 x 4 - 1300x 5 - (Q z - Q w ) = 0 

Similar to the mass balance, the energy balance constraints 
are also always defined as equality constraints. It may be 
noted that the energy balance is not normally carried out in 
LCA, mainly because of the difficulties in calculating the 
specific enthalpies for many of the substances and the 
enthalpy dissipated; the energy consumed in each system, 
is specified instead. 

The other constraints in this example include the limitation 
on the processing of the two categories of the feed, as de¬ 
fined by: 


Feed Category 1 availability: x, + x 2 S 1500 (A.3) 

Feed Category 2 availability: x 3 + x 4 + x 5 < 1200 (A.4) 

Therefore, the raw material availability can be defined as 
either equality or inequality constraints. The same is true 
for capacity, quality, heat etc. constraints. 


Finally, output of the functional units is subject to the fol¬ 
lowing constraints: 


Product 1 constraint: 

x 6 = 1000 

(A.5) 

Product 2 constraint: 

x 7 =1000 

(A.6) 

and 

x,,x 2 ,...,x 7 >0 

(A.7) 


The last constraint means that no activity can have negative 
values. Obviously, constraints (A.l)-(A.7) are equivalent to 
the general form of the LP model defined by eqns. (2) and 
(3) in the main text. The specific enthalpies in (A.2), for 
example, represent the factors of proportionality (or input/ 
output coefficients) a ci in eqn. (2). 

The next step in formulating a LP model is definition of the 
objective functions. Here, profit and environmental burdens 
are defined as the objectives on which the system will be 
optimised. Taking the feed and energy cost, and product value 
coefficients listed above, the profit objective in this example 
is defined as: 

F= Product value - feed cost - energy cost 
Profit objective: 

F = 150x 6 + 120x 7 -1 10x, — 120x, — 130x,- 

, ' (A.8) 

110x 4 — 1 15 x 5 — 3 - 10~ 3 (x 6 + x 7 ) 

Equation (A.8) is equivalent to eqn. (1) defined in the main 
text. 

The burden is defined as the sum of the individual contribu¬ 
tions from each activity from cradle to the point of their 
entry into the system. In addition, the system itself produces 
a certain amount of the burden, equal to 5 kg for the total 
production of 2000 kg of products. The individual contri¬ 
butions to the total burden, expressed per unit level of activ¬ 
ity, are listed below: 

X, X, Xj x„ X, x 6 +x. 

Burden (kg/kg) 0.01 0.015 0.01 0.03 0.02 0.0025 

Therefore: 

Burden objective: 

B = O.Olx, +0.015x, -t-O.Olx, + 0.03x 4 +0.02 x 3 + 0.0025(x 6 + x 7 ) (A.9) 
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This equation corresponds to the general form of the bur¬ 
den objective function defined by eqn. (4). Similarly, if the 
analysis were carried out at the Impact Assessment level, the 
burden obiective function would be replaced by the impact 
objective, as defined by eqn. (5). 

System optimisation 

As mentioned above, LP optimisation can be used to iden¬ 
tify both economic and environmental optima of the sys¬ 
tem. If the focus of the analysis is the former, then the sys¬ 
tem in this example is optimised on the profit objective de¬ 
fined by (A. 8), subject to constraints (A.l )-(A.7), to give the 
solution shown in Table A.l. Total profit of £50000 is 
achieved by using feeds x, and x 4 only, while the other feeds 
are not being used in the system. The active constraints for 
this operating state of the system are defined by eqns. (A.3) 
and (A.5)-(A.6). The constraint (A.4) is non-active and has 
a slack value, i.e. unused amount of feed, of 700 kg. 

The total burden is also calculated at the solution of the 
model and is equal to 35 kg. This is equivalent to a conven¬ 
tional way of calculating the burdens in the Inventory stage, 
as most economic systems are operated around the point of 
optimum profit. Moreover, the marginal values are also cal¬ 
culated with the total burden to give the marginal allocation 
coefficients, as discussed in the main text. As only active 
constraints have non-zero marginal values, this means that 
the burdens are allocated among the Product 1, Product 2 
and the Feed category 1 availability constraints. Their re¬ 
spective values (not shown in Table A. 2) arc 0.0325, 0.0325 
and -0.02. This means that an increase in the production of 
either of the two products by one kg would increase the 
total burden by 0.0325 kg. However, as the marginal value 
of the Feed category 1 availability constraint has negative 
value, an increase of 1 kg in the feed availability would cause 
a 0.02 kg decrease of the burden. Thus, this simple example 
demonstrates that the burden does not have to be allocated 
fully to the products, but to the other constraints as well, 
depending on what determines the system operation. 


Table A.l: Optimisation results for the LP example 


Optimised on 

Activity Profit Burden 


Profit, F (£) 

50000 

26000 

Burden, B (kg) 

35 

25 

Feed 1, x, (kg) 

1500 

800 

Feed 2, x 2 (kg) 

0 

0 

Feed 3, x 3 (kg) 

0 

1200 

Feed 4, x 4 (kg) 

500 

0 

Feed 5, x 5 (kg) 

0 

0 

Product 1. Xg (kg) 

1000 

1000 

Product 2, x 7 (kg) 

1000 

1000 

Heat, Q (MJ) 

200 

320 


If the goal of the study is to identify options for environ¬ 
mental improvements as part of Improvement Assessment, 
the system is optimised on the burden (or impact) objective 
function. The results of this optimisation are also shown in 
Table A.l. For the same output of the products, the system 
operation is now determined by the active constraints (A.5) 
and (A.6). The total burden is equal to 25 kg, which repre¬ 
sents an improvement of 29% over the burden obtained in 
the optimisation on the profit. This reduction is achieved by 
using 700 kg less of the feed x, compared to the previous 
case and using x, instead of x 4 . However, the profit is re¬ 
duced by about 48%. In addition, the total energy used in 
the system increases from 200 MJ to 320 MJ. Therefore, in 
this particular case, the reduction in one burden leads to an 
increase in another. To help resolve these conflicting situa¬ 
tion, the system can be optimised on a number of environ¬ 
mental and socio-economic objectives function. The result¬ 
ing Pareto or non-inferior surface can then be used to trade¬ 
off the environmental and economic performance between 
the points of minimum burden and maximum profit to iden¬ 
tify the BPEO not entailing excessive costs. 


Erratum 


Int. J. LCA (5) 266 - 272 (1998) "Einstein's Lessons for Energy Accounting in LCA" 
by Rolf Frischknecht, Reinout Heijungs and Patrick Hofstetter 
on page 286 in the headline of section 4: The correct text is "E = me 2 " 
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